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Abstract 

We perform the analysis of existing light-targets deep-inelastic- 
scattering (DIS) data in the leading-order (LO), next-to-leading-order 
(NLO), and next-to-next-to-leading-order (NNLO) QCD approxima¬ 
tions and extract PDFs simultaneously with the value of the strong 
coupling constant ag and the high-twist contribution to the structure 
functions. The main theoretical uncertainties and experimental uncer¬ 
tainties due to all sources of experimental errors in data are estimated, 
the latter generally dominate for the obtained PDFs. The uncertainty 
in Higgs boson production cross section due to errors in PDFs is ~ 2% 
for the LHC and varies from 2% to 10% for the Fermilab collider un¬ 
der variation of the Higgs boson mass from 100 GeV to 300 GeV. For 
the IF-boson production cross section the uncertainty is ~ 2% for the 
both colliders. The value of = 0.1143 ± 0.0014(exp.) is 

obtained, while the high-twist terms do not vanish up to the NNLO 
as required by comparison to data. 

PACS numbers: 13.60.Hb,06.20.Jr,12.38.Bx 
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1 Introduction 


The deep-inelastic-scattering (DIS) data are used for extraction of the par- 
ton distribution functions (PDFs) in nucleon starting from the advent of the 
Bjorken scaling phenomena (for the most recent studies see Refs. [|], ||, |]). 
However since new DIS data appear permanently and the theoretical ba¬ 
sis for their analysis is being developed alongside, those analysis should be 
updated in order to improve the experimental and theoretical errors on the 
extracted PDFs. In this paper we give such update of the next-to-leading- 
order (NLO) QCD analysis of Ref. aimed to extract PDFs, the high-twist 
(HT) contribution to the structure functions F 2 ,l, and the value of strong 
coupling constant Og from the existing data on DIS of charged leptons off the 
proton and deuterium targets. One of the improvements of our analysis is 
replacement of data by the HI [Q and the ZEUS 0 collaborations by their 
more recent data 1,0 obtained in the 1996-97 run of the HERA collider, 
when these experiments accumulated about 5 times more events than before. 
Another improvement is account of the NNLO QCD corrections, which have 
been known with a good precision after calculation of the Mellin moments of 
the NNLO splitting functions up to 14-th ^j. This new input signihcantly 
decreases the uncertainty of PDFs due to the higher-orders (HO) QCD cor¬ 
rections 1,0 and together with more precise data allows to produce the 
NLO and NNLO PDFs sets with reduced total uncertainty. 

The particular features of the analysis [Q as compared to the global fits 
of Refs. [O, O are: 


• data on the differential cross sections are used that allows to use ex¬ 
perimental constraints on the structure function iA; 

• corrections for the effects of target mass (TM) and for the nuclear 
effects in deuteron are applied to the data, the both were calculated in 
the hts iteratively using the current PDFs; 

• comprehensive account of the experimental errors and their correlations 
is performed; 

• impact of main sources of the theoretical errors is analyzed. 


All these features are also kept in the present analysis. 
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Table 1: The numbers of data points (NDP) for separate experiments and the corresponding renormalization 
shifts ^ obtained in different hts. 


proton deuterium 


Experiment NDP 

e[%] 

NDP 

e[%] 




LO 

NLO 

NNLO 


LO 

NLO 

NNLO 

SLAC-E-49A 

58 

2.9 ± 1.4 

1.6 ± 1.4 

1.1 ± 1.4 

58 

0.0 ±1.3 

-1.3 ± 1.3 

-1.7± 1.3 

SLAC-E-49B 

144 

3.2 ± 1.4 

1.8 ± 1.4 

1.2 ± 1.4 

135 

0.5 ± 1.4 

-0.8 ± 1.4 

-1.3 ± 1.4 

SLAC-E-87 

90 

2.9 ± 1.4 

1.8 ± 1.4 

1.2 ± 1.4 

90 

0.7 ± 1.3 

-0.5 ± 1.3 

-0.9 ± 1.3 

SLAC-E-89A 

66 

5.6 ± 1.9 

4.0 ± 1.9 

3.1 ± 1.9 

59 

2.3 ±2.0 

0.5 ±2.0 

-0.3 ±2.0 

SLAC-E-89B 

79 

2.6 ± 1.3 

1.2 ± 1.4 

0.7 ± 1.4 

62 

-0.2 ± 1.3 

-1.4 ± 1.3 

-1.9 ± 1.4 

SLAC-E-139 





16 

1.6 ± 1.3 

0.2 ± 1.3 

-0.2 ± 1.3 

SLAC-E-140 





26 




BCDMS 

351 




254 




NMC(90 GeV) 

46 

-1.5 ± 1.6 

-1.4 ±1.5 

-1.3 ± 1.6 

46 

-2.6 ± 1.5 

-3.0 ± 1.5 

-2.9 ± 1.5 

NMC(120 GeV) 

59 

-0.2 ± 1.5 

0.2 ± 1.5 

0.4 ± 1.5 

59 

-1.3 ±1.5 

-1.6 ± 1.4 

-1.5 ± 1.5 

NMG(200 GeV) 

66 

2.1 ± 1.5 

1.9 ± 1.4 

2.0 ± 1.4 

66 

0.0 ± 1.4 

-0.3 ± 1.4 

-0.3 ± 1.4 

NMG(280 GeV) 

74 

1.3 ± 1.4 

0.9 ± 1.4 

0.8 ± 1.4 

74 

-1.1 ± 1.4 

-1.4 ± 1.4 

-1.4 ± 1.4 


HI 135 

ZEUS 161 



2 Theoretical and experimental input of the 
fit 

The model for data description is based on perturbative QCD with phe¬ 
nomenological parameterization of the LT and HT contributions to the struc¬ 
ture functions The analysis was performed in the MS scheme with the 
number of flavors fixed at 3. The LT PDFs are parameterized at Qq = 9 GeV^ 
in the following form: 


xuy{x,Qo)= + 72 0 ;), 

(1) 

xusix,Qo) = - 0 ;)^'=+ 

iVs 

(2) 

xdyix, Qo) = -^X°-’^{1 - 

(3) 

xds{x,Qo) = 

(4) 

a;ss(o;,Qo) = 

(5) 

xG{x, Qo) = - o:)^°(l + 7 p\/x + 7 ^ 0 ;), 

(6) 

where u, d, s, G are the up, down, strange quarks, and gluons distributions 
respectively; indices V and S correspond to the valence and sea quarks. The 
parameters and Aq were calculated from the other parameters us¬ 

ing conservation of the partons momentum and the fermion number. The 
normalization parameter is selected in such way that Hg gives the to- 


tal momentum carried by the sea quarks. The parameter rjs was fixed at 
0.42, in agreement with the results of NuTeV collaboration |^. The rest 


of parameters coming to Eqns.(|I|-§D were fitted to the data, we convinced 
that no extra parameters are needed to improve the data description. The 
HT contributions to the fitted structure functions F 2 ^l were parameterized 
in additive form 

TT _ 77LT,TMC I H2;l{x) 

J^2,L — -C, T -h 


2,L 




4 





Table 2: The values of x^/NDP for different fits and of the averages of 
residuals R for the NLO £t together with their standard deviations AR for 
separate data sets used in the analysis. 


NDP xVNDP R{AR) 


Experiment 


LO 

NLO 

NNLO 


SLAG-E-49A 

115 

0.54 

0.52 

0.52 

-0.04(0.23) 

SLAG-E-49B 

279 

1.27 

1.20 

1.21 

0.23(0.29) 

SLAG-E-87 

180 

0.96 

0.92 

0.91 

0.01(0.37) 

SLAG-E-89A 

125 

1.35 

1.33 

1.36 

-0.18(0.45) 

SLAG-E-89B 

141 

0.94 

0.82 

0.81 

0.47(0.49) 

SLAG-E-139 

16 

0.81 

0.58 

0.62 

-0.07(0.43) 

SLAG-E-140 

26 

1.31 

0.94 

0.90 

0.17(0.86) 

BGDMS 

605 

1.15 

1.13 

1.12 

-0.06(0.68) 

NMG 

490 

1.35 

1.24 

1.24 

-0.06(0.39) 

HI (96-97) 

135 

1.39 

0.97 

1.09 

0.42(0.55) 

ZEUS(96-97) 

161 

1.54 

1.32 

1.28 

-0.55(0.64) 

TOTAL 

2274 

1.20 

1.11 

1.11 

0.01(0.22) 


where account of the TM correction^ and 

functions H{x) are parameterized in the piece-linear form and htted to the 
data. 

The code used to evolve these distributions was used earlier in the anal- 

m 


ysis of Ref. 


This code was checked against Les Houches benchmark of 
Ref. and demonstrated accuracy much better than experimental precision 
of the analyzed data. A brief description of the data sets used in the analysis 
is given in Table |l|. We cut data with Q'^ < 2.5 GeV^, > 300 GeV^, 
and X > 0.75 in order to minimize effects of the higher-order QGD correc¬ 
tions, electroweak interference, and uncertainties in the TM and deuteron 
corrections correspondingly. All experimental nncertainties in data released 
by authors were accounted for in the analysis, including statistical errors 
and correlated systematical errors (in particular the general normalization 
errors). 

The total normalization of the NMG and most of the SLAG data subsets 


^See Ref. 101 for detailed formula. 
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Figure 1: The experimental (statistical and systematical) errors bands for 
the PDFs obtain in the LO (dots), NLO (dashes), and NNLO (full) hts. 
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used in our analysis were not determined in the experiments, instead they 
were htted by authors to the results of other experiments. Since now we 
have more DIS data than it was used in those hts, we perform independent 
renormalization of the NMC and SLAG data subsets that allows for more 
precise determination of the corresponding normalization factors. Following 
the approaches adopted by those experimental groups we introduce an addi¬ 
tional normalization parameter for each beam energy and target in the NMC 
data set and for each experiment and target in the SLAG data sets. The data 
from SLAC-E-140 experiment, which was performed with particular atten¬ 
tion paid to the accurate determination of the absolute normalization, were 
not renormalized - we accounted for its general normalization error given by 
the authors in the same way as the other correlated systematic errors. The 
same approach was used for treatment of the general normalization error in 
the BCDMS, HI, and ZEUS data. The normalization factors for the SLAG 
and NMC data are given in Table |I| for the hts performed in the LO, NLO, 
and NNLO approximations. In the LO ht the shifts for the SLAC deuterium 
sets deviate off zero by about two standard deviations, while for the NNLO 
ht they vanish, i.e. account of the higher-order corrections decreases tension 
between the diherent data sets (evidently, this ehect is less pronounced for 
the NMC data than for the SLAC one since the later correspond to the lower 
values of Q, where these corrections are more important). In average, renor¬ 
malization of the SLAC data sets in the NLO ht is ~ +2% for the proton and 
~ —0.5% for the deuterium targets. For the NMC data the corresponding 
values also are diherent and, besides, vary with the beam energy that justihes 
application of the hexible renormalization scheme for this experiment. Note 
that in particular due to this hexible renormalization we obtain satisfactory 
description of the NMC data in the NLO: x^/NDP=1.24 in contrast with 
value 1.5 obtained in the ht of Ref. |n|. The values of y^/NDP obtained in 
diherent orders of QCD are given in Table |^. One can see that quality of the 
ht is approximately the same for the NLO and NNLO hts and much worse 
for the LO ht. 


3 PDFs and their experimental errors 

The PDFs obtained from the hts in diherent pQCD orders with their to¬ 
tal experimental uncertainties are given in Fig.|^ the PDFs parameters are 
given in Table One can see that the shift between the NLO and NNLO 
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Q^=9 GeV^ 



Figure 2: Impact of changes in the data set on the NLO PDFs and their 
errors. Full lines: relative experimental errors in our PDFs; dashes: the 
same for the PDFs of Ref. 0; dots: the same for the CTEQ6 PDFs. 



















PDFs is generally within their experimental errors excluding the sea quarks 
distribution at small x. The LO gluon distribution at small x and valence 
quarks distributions at moderate x are very different from the corresponding 
NLO and NNLO PDFs. This difference cannot be completely ascribed to 
the re-definition of PDFs in different orders since data are poorly described 
by pQCD in the LO. Fitted LO PDFs are evidently distorted due to inac¬ 
count of the HO corrections and this distortion should be taken into account 
in calculation of other hard processes in the LO approximation. The PDFs 
errors consecutively fall from LO to NNLO due to additional terms of the 
perturbative expansions put additional constraints on the fitted values. This 
effect is most pronounced for the gluon distribution at large x. The com¬ 
parison with experimental uncertainties in the NLO PDFs obtained from 
our earlier analysis of Ref. is given in Fig.^. The main improvement is 
in precision of the gluon distribution at small x due to more precise HERA 
data have been used. The total experimental errors can be conventionally 
split into 3 components: the uncorrelated, general normalization, and cor¬ 
related errors (excluding normalization ones). The normalization errors can 
also be considered as correlated ones, but we separate them for the purpose 
of comparison since very often the normalization errors are accounted for in 
a peculiar way or not accounted for at all. In our analysis the normalization 
errors for BCDMS, HI, and ZEUS data sets are included in the general co- 
variance matrix, while for SLAG and NMC data sets they appear through 
the error matrix of the fitted parameters. Ratios of the correlated errors to 
the uncorrelated ones and of the normalization errors to the correlated ones 
for different PDFs are given in Fig.^. One can see that for most of PDFs the 
correlated errors are more important at small x. The normalization errors 
are largest at large x; in this region they give main contribution to the total 
experimental uncertainties in PDFs. 
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Q^=9 GeV^ 



Figure 3: Ratios of the correlated systematic errors to the uncorrelated ones 
(upper panel) and of the normalization errors to the correlated ones (lower 
panel) for different PDFs (full curves: gluons; dashes: sea quarks; dots: d- 
quarks, dashed-dots: u-quarks). 
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Table 3: The values of the htted parameters of the PDFs. 




LO 

NLO 

NNLO 





(A±B)/2 

B 

Valence 


0.551 ±0.085 

0.700 ±0.030 

0.725 ±0.027 

0.725 ±0.026 


K 

3.672 ± 0.042 

3.920 ±0.050 

4.024 ±0.050 

4.011 ±0.050 


72 

3.0 ±1.7 

1.14 ±0.39 

1.05 ±0.35 

1.01 ±0.34 



0.639 ±0.047 

0.722 ±0.082 

0.772 ±0.074 

0.818 ±0.062 


bd 

4.48 ±0.23 

4.94 ±0.12 

5.14 ±0.15 

5.22 ±0.18 

Glue 







O-G 

-0.302 ±0.021 

-0.146 ±0.018 

-0.128 ±0.021 

-0.096 ±0.020 


be 

5.3 ±5.1 

8.2 ± 1.3 

9.4 ±1.1 

9.8 ± 1.0 


7 ? 

-1.94 ±0.83 

-3.76 ±0.45 

-3.84 ±0.52 

-4.18 ±0.45 


7 ^ 

2.8 ±5.3 

7.7 ± 1.7 

8.6 ±1.7 

9.2 ± 1.5 

Sea 

71s 

0.161 ±0.013 

0.160 ±0.011 

0.1591 ±0.0090 

0.1545 ±0.0071 


^sd 

-0.1980 ±0.0057 

-0.1968 ±0.0048 

-0.2092 ± 0.0044 

-0.2086 ±0.0043 


bsd 

3.8 ± 1.1 

5.1 ± 1.4 

5.6 ±1.3 

6.5 ± 1.1 


Vn 

1.46 ±0.15 

1.16 ±0.12 

1.12 ±0.11 

1.036 ±0.099 


bsn 

9.2 ± 1.8 

10.16 ±0.93 

10.49 ±0.90 

10.43 ±0.89 



Impact of each data set on the precision of our PDFs is demonstrated in 
Fig.^. The HERA data are crucial for determination of the gluon and quarks 
distributions at small x; the BCDMS data constraint the quark distributions 
at large x; the NMC data are essential for determination of the anti-quarks 
and d-quark distributions at large x; and the SLAG data improve precision 
of the d-quark and gluon distributions at large x. In summary none of the 
subsets can be dropped out without deterioration of the PDFs precision. We 
also checked how sensitive are the errors in PDFs to the inclusion in £t the 
ZEUS data of Ref. [0 with Q‘^ > 250 GeV^, the HI data of Ref. 

> 150 GeV^, and the FNAL-E-665 data of Ref. 
these data our regular kinematical cuts 


|l6| with 


rfi 


If one applies to 
> 2.5 GeV^ and x < 0.75, the 


relative suppression of the PDFs errors after each of these data sets has been 
included does not exceed 0.9 and therefore they are useless for the PDFs 
determination. 

The novel HERA data used in our £t is confronted with the htted curves 
in Figs.|,|. For comparison we show in the same hgures the calculations 
based on our earlier PDFs of Ref. ^j. The experimental error bands for 
both sets of PDFs overlap that proves stability of the analysis. If the £t 
is unbiased, the fitted curve and data should coincide modulo statistical 
fluctuations. These fluctuations are quantihed by the average of residuals 
R =< {f — y)/e >, where / is the fitted function, y is measurement, and 
e is the total experimental error in data. The standard deviation of this 
average AR is ~ 1/^/NDP for data set with uncorrelated points. If data 
are correlated, coherent shift of the htted curve within the value of common 
systematic error in data is allowed and value of AR gets larger, up to 1 for 
the case of strong correlations. Therefore the value of AR can be used as 
a crude indicator of signihcance of correlations in data set. In our analysis 
AR is largest for the BGDMS and SLAG-E-140 data (see Table ^). Since 
magnitudes of R are well within AR both for all subsets and for the total 
set we conclude that our fit is unbiased. Despite the systematic errors are 
not necessarily Gaussian distributed the value of x/^DP = 1.1 obtained in 
our £t proves that this is good approximation in our case. Distribution of 
the diagonalized residuals (DRs), which are the components of the vector of 
residuals {f — y)/e multiplied by the square root of the inverse covariance 
matrix is plotted in Fig.^ If the errors obey the Gaussian distribution the 
distribution of DRs have to be normal |T^ and as one can see in Fig.^this is 
the case for our £t. Therefore all errors are Gaussian in good approximation 
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Q^=9 GeV^ 



Figure 4: Relative increase of the experimental errors in our NLO PDFs due 
to rejection of different data sets from the fit (full curves: BCDMS data are 
rejected; dashes: HERA ones; dots: SLAG ones; dashed-dots: NMC ones). 
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Figure 5: The experimental error bands for calculated using our NLO 
PDFs (full curves) compared to the HI data used in the fit. Figures at the 
curves are values of in units of GeV^. For better view the factor of 
and vertical shifts are applied to the data points and curves. The same bands 
calculated using the A99 PDFs are also given for comparison (dashed curves). 
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Figure 6: The same as Fig.^ for the ZEUS data used in the £t. 
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Figure 7: Distribution of the diagonalized residuals in our NLO fit. The 
curve superimposed gives normal distribution. 

and consequently the errors in htted parameters are Gaussian tooQ. 


4 The PDFs theoretical uncertainties 

With experimental errors in PDFs permanently reduced due to to more rele¬ 
vant data appear the theoretical errors become more dominating contribution 
to the total PDFs uncertainty. In study of the importance of theoretical un¬ 
certainties the experimental errors serve as a natural scale. The variations 
below this scale can be neglected and, moreover, they just can occur due to 
statistical fluctuations in data. For this reason we analyzed only the the¬ 
oretical errors, which give effect larger or comparable to the experimental 
ones. 

^The width of the DRs distribution is 1.05 that is marginally larger than the width of 
the normal one; for this reason rescaling of the errors in spirit of the PDG’s procedure 
would lead to negligible effect. 
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Q^ = 0 Gev^ 



X 


Figure 8: Perturbative stability of the gluon distribution at small x. The 
corresponding RS error in the NLO and NNLO is given for comparison. 

4.1 Uncertainty due to higher order QCD corrections 

The uncertainty in PDFs due to the HO QCD corrections is not an un¬ 
certainty in full sense because the PDFs differ from order to order just by 
definition. Nevertheless the variation of PDFs with the pQCD order is of 
practical importance for estimation of the HO QCD uncertainties in the 
hard cross sections, which are calculated using the PDFs, despite the varia¬ 
tion of the cross section value due to the HO corrections to PDFs can be in 
some cases partially compensated by the HO corrections to the correspond¬ 
ing coefficient functions. A common approach used for estimation of the 
HO errors in the analysis of DIS data is to vary the QCD renormalization 
scale (RS) and factorization scale (FS) from to 4Q^. Evidently, since the 
range of this variation is arbitrary and no a:-dependence of the scale factors 
is assumed, this approach gives only semi-qualitative estimate on the HO 
uncertainty. Having the PDFs at different pQCD orders we can verify this 
approach comparing the real variation of PDFs from order to order with the 
HO errors estimated using the variation of the pQCD scales. Example of 
such comparison is given in Fig.H. One can see that the real change of the 
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gluon distribution from the NLO to the NNLO is much smaller than the RS 
error at small x. With the FS errors taken into account the difference would 
get even larger. Adjustment of the range of the scale variation would also 
not help since the signs of the NLO-NNLO difference and the RS error are 
opposite at X ~ 0.01. Note that the RS error in gluon distribution almost 
does not decrease from the NLO to the NNLO. This means that if the per¬ 
turbative expansion does converge the discrepancy between the NNLO RS 
error and the change of gluon distribution from the NNLO to the NNNLO 
is even larger. For this reason we recommend to use the difference between 
NLO and NNLO PDFs as estimate of the HO uncertainty motivated by the 
properties of the asymptotic expansions. The NNLO PDFs have an addi¬ 
tional source of errors due to possible variation of the NNLO kernel. The 
magnitude of these errors is largest at small x, but even in this region they 
are smaller than the experimental errors (see Fig.^. 


4.2 Uncertainties due to the heavy quarks contribu¬ 
tion 

4.2.1 Variation of the factorization scheme 


In our analysis contribution of the c- and 6-quarks was calculated in the hxed- 
flavor-number (FFN) scheme using the photon-gluon fusion mechanism 


with account of the NLO QCD correction [^. Since the heavy-quark con¬ 


tribution rises with Q rapidly the FFN scheme without resummation of large 
logs is irrelevant at small x and Q much larger than the quark mass . In 
order to overcome this difficulty the variable-flavor-number (VFN) scheme 
was proposed |^. In the VFN scheme the heavy quarks are considered as 


massless and are included in the QCD evolution at Q = mc,b,t, where mc,b,t 
are the masses of c-, b-, and t-quark correspondingly. The VFN scheme does 
include resummation of large logs, but at the same time evidently overesti¬ 
mates the heavy-quarks distributions at Q close to the thresholds. Never¬ 
theless this scheme is widely used in practice due to it allows for to greatly 
simplify the calculation of processes with the heavy quarks involved. The 
FFN scheme is appropriate in our case since the heavy-quarks contribution 
is concentrated at small x, while the small x data used in our analysis corre¬ 
spond to the values of Q, which are not very far from the threshold of heavy 
quarks production. Nevertheless in order to lighten implementation of our 
PDFs in different calculations we generated the VFN PDFs using our FFN 
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Q'=9 GeV' 



10 10 ^ 10 ^ X 


Figure 9: The NNLO PDFs obtained with the different choices of the NNLO 
kernel (full lines: la experimental bands for the nominal fit; dashes: fit with 
the NNLO kernel modihed within the allowed uncertainty given in Ref.[0). 


19 




xc 


xb 




Figure 10: The c- and 6-quarks distributions in the VFN scheme (dashed 
curves) compared to the corresponding FFN contributions to the structure 
function F 2 weighted with the quarks charges (full curves). The hgures at 
the curves give the values of in GeV^; the range of x for a given is 
typical for the processes studied at the Fermilab and LHC colliders. 


PDFs as input. For this purpose we evolved our PDFs from our boundary 
value of Qq = 9 GeV^ to Q = rric in the FFN scheme and used the obtained 
PDFs as the boundary condition for the VFN evolution. The obtained VFN 
distributions of c- and 6-quarks in the NLO are given in Fig.^ in comparison 
with the corresponding contributions to the structure function F 2 calculated 
in the FFN scheme and weighted for the quarks charges. At large Q the FFN 
scheme gives results similar to the VFN ones at the values of x feasible by 
future hadron colliders. At low Q the VFN scheme gives sizeable overesti¬ 
mation of the heavy-quark distributions, especially for the 6-quarks. For this 
reason the certain caution is necessary in use of the VFN PDFs at small Q 
in this connection). 


(see also Ref. [22 
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4.2.2 Uncertainty due the error in heavy quarks masses 

This uncertainty is estimated as the variation of PDFs under the change of 
rric from 1.5 GeV to 1.75 GeV, i.e. within the error in its value estimated 
from other measurements. The magnitude of corresponding uncertainties in 
PDFs decreases from the LO to the NNLO since the HO PDFs are smoother 
at small x and variation of the photon-gluon fusion cross section gets less 
sensitive to the shift of the kinematical variables. The ratio of its magnitude 
to the PDFs experimental errors in the NLO is given in Fig.^. The ratio is 
generally less than 1.5 being most important for the sea quarks distribution 
at small x. The similar ratio for uncertainty due to variation of the 6-quark 
mass is negligible for all PDFs. 

4.3 Uncertainty due to variation of the strange sea pa¬ 
rameterization 

The s-quark distribution equally contributes to all structure functions of the 
charged leptons DIS and cannot be determined from these data due to lack 
of constraints. Most precise information about this distribution comes from 
the neutrino nucleon DIS data. Its shape was determined by the NuTeV 
collaboration to be comparable to the average of the u- and d-quarks distri¬ 
butions and the total suppression factor of the s-quark distribution rjs was 
found to be 0.42 ± 0.07(stat.) ±0.06(syst.) We estimate the uncertainty 
in PDFs due to error in the s-quark distribution as change of the PDFs under 
variation of r]s from 0.42 to 0.52. This uncertainty marginally changes from 
the LO to the NNLO and does not exceed the PDFs experimental errors (see 

Fig.|n]). 

4.4 Other theoretical errors 

The list of possible theoretical errors in PDFs is rather conventional. For 
example the errors due to HT contribution can be estimated as shift of the 
fitted PDFs under addition to the power-like terms with the x-shape moti¬ 
vated by a certain model considerations to the htted function. Alternatively, 
it can be accounted for by htting possible HT contribution in the model- 
independent form; in this case the error due to the HT is combined with the 
total experimental error in PDFs. Since the later approach is less subjective 
we use it in our analysis. Similarly, we do not treat the error due to the 
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Q^=9 GeV^ 



Figure 11: The ratio of magnitude of the uncertainty due to errors in the 
c-quark mass (upper panel) and the strange sea suppression factor (lower 
panel) to the experimental errors in NLO PDFs (full lines: gluon distribution; 
dashes: non-strange sea one; dots: d-quark one; dashed-dots: u-quark one). 
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uncertainty in value of Ofs as a separate theoretical error since we fit si¬ 
multaneously with other parameters. In our previous analysis of Ref. 0 we 
considered the error due to uncertainty in the heavy-quarks thresholds in the 
QCD renormalization group equation for as- Now we omit this error because 
it is generally negligible as compared to the experimental errors and other 
theoretical errors. The errors due to uncertainty in the nuclear effects in 
deuterium, which also was considered in Ref. Q, will be analyzed in details 
elsewhere 
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5 High-twist contribution and the value of a. 


In the analysis of DIS data with relatively low Q one needs to account for 
the HT terms due to these terms are essential part of the operator product 
expansion ||2^. The HT terms cannot be accurately calculated from the 


existing theory of strong interaction and we have to account for them in 
our analysis phenomenologically. From comparison with the data we observe 
non-negligible contribution from the HT terms since account of such terms 
improves quality of the £t (e.g. in the NNLO fit the value of rises by 
330, if we £x the HT terms at 0). At the same time the observed HT terms 
can come just from the possible defects of the htted model, in particular 
due to inaccount of the HO QCD corrections. Indeed, with the inclusion of 
the NLO correction in the analysis of data on the uN DIS the observed HT 
terms decrease and eventually vanish in the NNLO [^. We also observe 
decrease of the HT terms with the pQCD order, but they do not vanish 
completely in our case (see Fig.|^ and Table^. Extrapolation of the LO- 
NLO-NNLO results for the HT terms allows to expect that they should not 
vanish in the NNNLO also. This is in line with the conclusion of Ref. 
that the HO corrections can successfully simulate power behavior in narrow 
region of x close to 0.1 only. This is also in agreement with the results of 
Ref. (the effects of soft gluon resummation considered in Ref. are 
important at x > 0.8 only and cannot further decrease the value of HT terms 
in our case). We see no contradiction between our results and the results of 


Ref. 1^ in view of large errors in the HT terms extracted from the neutrino 
data, but do not support the conclusion about vanishing of the HT terms 
in the NNLO. The obtained HT contribution to F 2 is maximal at a: ~ 0.6 
and at = 5 GeV^ is 10% of the LT term. This is much larger than the 
value of the relative experimental error in the F 2 measured in this region and 
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Table 4: Fitted coefficients of the HT contributions to the structure functions 
F 2 ,l in the NNLO 


X 


Hdix) 

0. 

-0.034 ±0.022 

0.234 ±0.037 

0.1 

-0.056 ±0.017 

0.042 ±0.021 

0.2 

-0.019 ±0.011 

0.055 ±0.021 

0.3 

-0.003 ±0.010 

0.072 ±0.018 

0.4 

0.0077 ±0.0098 

0.002 ±0.017 

0.5 

0.0299 ±0.0082 

0.0138 ±0.0087 

0.6 

0.0469 ±0.0077 

0.0060 ±0.0091 

0.7 

0.0447 ±0.0062 

0.0233 ±0.0063 

0.8 

0.0398 ±0.0065 

0.026 ±0.012 


Table 5: Values of as(Mz) obtained in different orders of pQCD. 

LO 0.1301 ±0.0026(exp) 

NLO 0.1171 ±0.0015(exp) 

NNLO 0.1143 ±0.0014(exp) 


hence the account of the HT terms is indeed important. 

Values of as{Mz) obtained in different orders of pQCD are given in Ta¬ 
ble ^ For comparison in the NNLO analysis of the existing data for DIS of 
charged leptons off proton the value of as(Mz) = 0.1166 ± 0.0009 (exp) was 
obtained [^]. This value is comparable with our result, but one is to note 
that, despite the low Q data were used in the £t, the analysis of Ref. 
was performed with the HT terms hxed at zero. Meanwhile the values of as 
and H 2 at large x are strongly anti-correlated. From the trial £t with the 
constraint 772 ,l = 0 we obtain the value of ag(Mz) = 0.1215 ± 0.0003 (exp), 
i.e. their value is much larger and their error is much smaller than in the 
£t with the HT terms htted. From this comparison we conclude that in fact 
our result for ag is in disagreement with the result of Ref. 
value of ttg obtained in the analysis of Ref. 

NNLO, in contrast with our results. 


28[|. In addition. 


28i rises from the NLO to the 
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Figure 12: The HT contributions to the proton structure functions F 2 (upper 
panel) and the nucleon structure function Fl (lower panel) obtained in the 
different orders of pQCD. For better view points are shifted to left and right 
along the x-axis. 
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6 Comparison with other parameterizations 


Our NLO PDFs evolved in the VFN scheme are compared to the recent 
CTEQ and MRST parameterizations in Fig. |l^. The comparison is made 
at = 100 GeV^, which is the lowest scale for the hard processes at the 
hadron colliders. At this scale we observe the largest deviations of our PDFs 
from the CTEQ gluon distribution at large x, the MRST gluon distribution 
at small x, and the quark distribution of both groups at small x. The ob¬ 
served differences in particular can be due to the data on Drell-Yan process 
and jet production were used in the CTEQ and MRST analysis. To clarify 
this point thorough analysis of universality of the PDFs in different processes 
is necessary, but any way these differences can have impact on the calculated 
cross sections of hard processes at the hadron collider, e.g. on the K-factor 
for the Higgs boson production. The discrepancy with the MRST gluon dis¬ 
tribution, which is the lowest of the three considered PDFs sets, is especially 
important since at small x it gets negative at small Q. At larger Q the dif¬ 
ference between the different PDFs sets is smaller due to focusing properties 
of the QCD evolution [p^] . 


Recently the CTEQ collaboration estimated the errors in their PDFs and 
we can compare those estimates with the errors in our PDFs. The comparison 
cannot be completely sensible since due to technical difficulties CTEQ did 
not take into account the impact of the normalization errors in data. Indeed, 
the account of normalization errors essentially increase the point-to-point 
correlations for many subsets and leads to amplihcation of the numerical in¬ 
stability of the calculations. The value of AR, which is indicator of these cor¬ 
relations (see Sec.§), rises essentially if normalization errors are included into 
the covariance matrix. In particular, for the BCDMS data AR is 0.68/0.29 
with/without normalization errors taken into account. Meanwhile account 
of the normalization errors leads to increase of the PDFs uncertainties up to 
factor of 5 (see Fig.|j). In the analysis of CTEQ the value of x^/^DP for 
some data subsets deviates off 1 much more than this is statistically allowed. 
These fluctuations may be caused by several reasons, including inaccuracies 
in the model of data, such as inaccount of the HT terms, the target mass 
corrections, the nuclear corrections, and the HO QCD corrections. To avoid 
this problem CTEQ scaled their PDFs errors by the “tolerance” factor, which 
was selected equal to 10 in order to provide comparability of all data subsets. 
Effectively, the application of such tolerance factor leads to the crude account 
of all sources of uncertainties, including theoretical ones, since the scaled er- 
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Q"= 100 GeV" 
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Figure 13: One-standard-deviation bands for our NLO PDFs evolved in the 
VFN scheme (full lines) compared to the MRST2001 (dashes) and CTEQ6 
(dotes) PDFs. 
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Figure 14: One-standard-deviation bands for the Q-dependence of the mo¬ 
mentum carried by different parton species [q reads sum of the anti-quark 
distributions). 


rors cover the total range of possible variation of the htted parameters. In 
our analysis we need not introduce the tolerance factor since deviations of 
X^/NDP for separate data subsets off 1 £t to the possible statistical fluc¬ 
tuations (see Table H). The errors in CTEQ6 PDFs are compared to ours 
in Fig. 1^. For all PDFs, excluding the sea quarks distribution at a: ~ 0.3 
the CTEQ PDFs errors are larger than ours. The difference is especially 
pronounced for the gluon and sea distributions at small x and for the u- 
and d-quark distributions at large x. From this comparison follows that use 
of the variety of data for the PDFs extraction not necessarily leads to the 
decrease of the PDFs errors since improvement in experimental accuracy can 
be accompanied by the growth of theoretical errors (at least one is to take 
care about accurate estimation of the later). 
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7 Access to the PDFs 


The PDFs obtained in onr analysis are available throngh the WWW pageQ. 
We placed in this page the PDFs generated in the LO, NLO, and NNLO; 
in the VFN and FFN schemes, inclnding special sets with shifted values 
of the c-quark mass and the sea suppression factor to allow account of the 
corresponding theoretical uncertainties in PDFs. The PDFs are supplied by 
the experimental uncertainties, which can be conveniently taken into account 
in the calculation. This can be performed both using propagation error 
formula based on the PDFs derivatives and by generating random PDFs. 
The later way is most convenient in the Monte-Carlo calculations and was 
adopted in Les Houches interface for the PDFs with uncertainties [H. The 


Icr bands for the hrst Mellin moments of the parton distributions at different 
calculated using this code are given in Fig.|T^ The typical uncertainties in 
the momentum carried by different species is generally less than 0.01, being 
the largest for the gluon distribution. The experimental errors in typical 
parton luminosities L for the antiproton-proton collisions at ^/s = 2 TeV 
(Fermilab collider) and the proton-proton collisions at ^/s = 14 TeV (LHC 
collider) are given in Fig.|T^. In the kinematical range relevant for the studies 
of hard processes the relative errors in parton luminosities for the LHC are 
in the range of 1 -i- 20%; for the Fermilab collider they are comparable with 
the LHC ones at M < 100 GeV and are much larger at M > 100 GeV. 
The uncertainty in VF-boson production cross section due to errors in PDFs 
2% for the both colliders; for the Higgs boson production cross section 


IS 


it is ~ 2% for the LHG and varies from 2% to 10% for the Fermilab collider 
under variation of the Higgs boson mass from 100 GeV to 300 GeV. 


8 Conclusion 

The PDFs extracted from the DIS data only are precise enough: The esti¬ 
mated uncertainties in the parton luminosities including all sources of ex¬ 
perimental errors in data are < 20% in the kinematical region feasible at 
the LHG collider and < 60% at the Fermilab collider. The theoretical errors 
in PDFs are under control and generally do not exceed experimental ones, 
therefore the DIS PDFs can be used for conclusive searches of deviation from 

^ ht tp: / / sirius .ihep. su/'alekhin/pdfa02 / 
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Figure 15: The dependence of parton-parton luminosities for the LHC 
(full curves) and the Fermilab collider (dashes) on the produced mass M 
(qq : Luu + Ldd + Ldu; qq : -f^udT-^da; (q + q)G- L^^Q-\-L^q-\-L^g L^g)■ 
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the Standard Model in other processes and/or checks of the nniversality of 
PDFs. 
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